Measuring the quantum states of a mesoscopic SQUID using a small Josephson 

junction 
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We have experimentally studied the energy levels of a mesoscopic SQUID using inelastic Cooper- 
pair tunneling. The tunneling in a small Josephson junction depends strongly on its electromagnetic 
environment. We use this fact to do energy level spectroscopy of a SQUID-loop by coupling it to a 
small junction. Our samples with strong quasiparticle dissipation are well described by a model of a 
particle localized in one of the dips in a cosine-potential, while in the samples with weak dissipation 
we can see formation of energy bands. 
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The Josephson junction, despite its simple structure 
has proven to be surprisingly versatile and new appli- 
cations are found in quantum computing and nanoelec- 
tronics Q, @- The devices are based on the quantum 
mechanical behavior of the superconducting phase vari- 
able , which has been previously studied with either rf- 
irradiation Q or during rapid current ramping 0, 0] . We 
are using a different probe, namely, an additional meso- 
scopic Josephson junction. Our scheme is based on the 
theory of phase fluctuations 0, 0], according to which 
Coulomb blockade in a single superconducting tunnel 
junction is strongly affected by its environment. Non- 
coherent Cooper pair tunneling is allowed only if energy 
is exchanged with the surroundings. Thus, this inelastic 
Cooper pair tunneling provides a good tool for observ- 
ing all kinds of environmental modes in a rather simple 
fashion 0. 

In this Letter we present detailed spectroscopic investi- 
gations on small SQUID loops, which are driven from the 
nearly classical limit (Ej/Ec ^> 1) deep into the quan- 
tum regime (Ej/Ec ~ 1). Our results yield evidence for 
higher energy bands of the macroscopic phase variable 
in a regime {Ej I Eg ~>. 1) where they have not been in- 
vestigated before [lfj. In addition, our experiment pro- 
vides the first verification that multiphoton transitions 
between the levels of a quantum mechanical harmonic 
oscillator play a prominent role in electron tunneling in 
a mesoscopic tunnel junction. 

As an energy detector in our measurement we use a 
voltage biased, superconducting tunnel junction which 
has a smaller size and critical current than the junction 
we want to study. For large (conventional) Josephson 
junctions the supercurrent is given by J = / c sin(</?), 
where I c is the critical current, which is related to the 
Josephson coupling energy Ej = HIc/(2e). The phase 
(p(t) = j 1 ^ jfV(t')dt' is defined as an integral of the 
voltage V across the tunnel barrier. For small junctions, 
where the charging energy Ec = e 2 /(2C) ^> Ej Cooper 



pair tunneling is inelastic and given by 



I(V) 



^ [P(2eV) 



P{-2eV)] 



(1) 



where P(E) is a function describing the probability of 
energy exchange between a tunnel junction and its elec- 
tromagnetic environment and depends on the impedance 
seen by the junction 

At low temperatures, the junction environment, i.e. 
the heat bath, is in its ground state and P{E) ~ for 
E < 0. Thus, the latter term in Eq. can be neglected 
and I(V) becomes directly proportional to P(2eV). The 
theory is valid for linear impedances constructed from 
lumped elements. Nevertheless, we argue that the idea 
of energy exchange can be generalized so that a discrete 
spectrum of energy levels in the environment will cause 
a set of discrete peaks in the /U-curve. Hence, the small 
detector junction can be used for spectroscopy. 

A Josephson junction can be described by the 
Schrodinger equation 
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+ -f-cos<^+ — ip )ip(tp)=0, (2) 
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where / is the current flowing through the junction. The 
current in our measurements always satisfies J <C Ic, so 
the tilt in the potential is negligible and setting 1 = 
in Eq. [2] leads to the familiar Mathieu-equation. The 
single junction Hamiltonian can also be used to describe 
a SQUID-loop, where the loop size is so small that the 
geometric inductance can be neglected and the loop is 
perfectly symmetric. The only difference is that Ej then 
depends periodically on an externally applied magnetic 
flux $ according to Ej = 2Ej" l9le \ cos(tt$/$ )| 2 , where 
$0 = V (2e) and E S j ngle is the Josephson coupling energy 
for a single junction. For large Ej/Ec, the particle is 
trapped in one of the potential wells. In this case, for 
currents I <^ Ic, the Josephson junction can be described 
by an inductance L = <j>o/(27r/ c ). Combined with the 
capacitance of the tunnel junction, the junction forms an 
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LC-oscillator with a characteristic resonance frequency of 



cj p = L/yLC = \/8EjE c /h. Consequently, a Josephson 
junction behaves like a harmonic oscillator with a level 
spacing of ui p . When Ej/Eq becomes smaller, the energy 
levels are not harmonic but they will depend on the shape 
of the cosine-potential. 

Depending on the environmental resistance seen by 
the Josephson junction, i.e. in our case "the environ- 
ment of the environment" , the junction can become 
completely delocalized and the whole periodicity of the 
cosine-potential has to be accounted for [ill fl2| . The 
eigenstates are then given by Bloch-functions ^n{f) = 
u n (if)e l(pq ^ 2e \ where q is the quasi-charge, n the band 
index and u n (tp) is a 27r-periodic function. This phase 
transition from the localized to delocalized state hap- 
pens when R > Rq, where Rq = /i/(4e 2 ), or 6.45 kfi 
[lfll Il4j | . In our measurement we need a clear voltage- 
bias and thus we have not fabricated any resistor close 
to the junction. The source of dissipation is, therefore, 
given by the quasiparticle resistance of the probe junc- 
tion. This changes the periodicity of the wave-functions 
from 2ir to Air and each band is split into two [l^ . 

We have carried out experiments with different circuit 
configurations; both 2- and 4- lead measurements includ- 
ing 1, 2 or 4 SQUID (s) coupled to a small detector junc- 
tion. We will here describe measurements of two differ- 
ent samples; a 4-SQUIDs sample, with four leads and a 
1-SQUID sample with just two leads. A scanning elec- 
tron micrograph of the 4-SQUID sample, together with a 
schematic drawing of it, is shown in Fig. 2] The SQUID 
configuration allows us to change the energy levels of the 
measured system and enables us to resolve the resonances 
due to the SQUID (s) from other resonances in the envi- 
ronment. The critical current, or equally, the value of Ej 
could be tuned to less than 1% of the maximum, which 
shows that our SQUIDs were very homogeneous. The 
samples were made from aluminum with e-beam lithog- 
raphy and 2-angle evaporation in an UHV chamber. 




The four wire setup facilitates the determination of cir- 
cuit parameters. The important parameters are Ej and 
Ec, or rather their ratio. The Ambegaokar-Baratoff (A- 
B) formula, Ej = Htt A/ (4e 2 Rt) , was used to find Ej 
from the normal state resistances, Rt, while the capaci- 
tances where estimated from the junction areas (see Fig. 
1) using a value of 45 fF/^m 2 15]. The BCS-gap, A, 
was about 215 fieV in our samples. The experimental 
parameters for the different circuits are summarized in 
Table 1. 

The samples were mounted into an rf-tight copper en- 
closure and cooled down to 80 mK with a plastic dilu- 
tion refrigerator. The measurement leads were filtered 
using 0.7 m long sections of Thermocoax. Mini-Circuits 
rf-filters with a cut-off frequency of 1.9 MHz were em- 
ployed on the top of the cryostat, at room temperature. 

Fig. |21 displays the measured TV-curve for zero mag- 
netic flux, or maximum Ej, for the 4-SQUID sample 
together with an IV-curve simulated with P(i5)-theory. 
The locations of the peaks were found to depend only 
on the magnetic flux, not at all on the gate voltages. 
To properly identify the energy levels of the SQUID we 
measure IV-curves for different magnetic fields. The peak 
positions as function of applied flux are shown in Fig|3| 

The width of the resonance peaks (about 4 /zeV) is 
smaller than ksT = 7fieV. The width is therefore, cither 
intrinsic or given by external noise. Our peak widths are 
thus comparable to or even smaller than what has been 
observed in similar spectroscopic studies [lfij ] . 

The peak structure can be qualitatively explained with 
a 3-resonator model, where one resonator represents all 
the SQUIDs and the two other come from the rest of the 
measurement circuitry; bonding wires and pads. The pa- 
rameters of the simulation and simulation circuitry are 
found in Fig. [5J The parameters for the SQUID were 
taken from independent measurements but the parame- 
ters of the two other resonator circuits were fitted to IV- 
curve. The resistances used in the simulation represent 
the broadening of peaks due to dissipation and noise. 

The P(_E)-function in Eq. was calculated using the 
integral equation approach presented in Ref. ITtI The 
comparison of the IV-curve with the simulation clearly 
shows that multiphoton excitations are present in the 
experiment. However, the energy levels are not equally 
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FIG. 1: SEM micrograph of a sample with 4 SQUIDs. The 
probe junction has an area of 100 x 100 nm 2 and the SQUID 
junctions 150 x 550 nm 2 In the samples covered in this paper, 
additional gates leads were available for the islands. Inset 
shows the schematic of the circuit in the 4- lead measurement. 



TABLE I: Parameters for the 4-squid and 1-squid samples. 
Energies are given in units of fieV. The Ambegaokar-Baratoff 
values for Ej are given in parentheses. 



Sample 


Rt (Ml) 


C(fF) 


Ej 


E C 


Ej/Ec 


4-SQUID (detector) 


166 


0.5 


3.6 


160 


0.023 


4-SQUID (squid) 


2.5 


7.6 


544 (272) 


10.5 


51.8 


1-SQUID (detector) 


70 


0.8 


8.5 


100 


0.08 


1-SQUID (squid) 


3.5 


5.7 


422 (188) 


14 


30.1 
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FIG. 2: IV-curve for the 4-SQUID sample at maximum Ej 
and the circuit model used in simulation. The full line shows 
the experimental curve while the shaded area shows the sim- 
ulated curve. The different excitations are denoted as (n,k,l), 
where n,k and 1 are the number of quanta excited. The first 
index gives the resonance due to SQUIDs and the two other 
indices are due to other resonances in the circuit. The pa- 
rameters used in the simulation are Criet = 0.5 fF, Ci= 4 fF, 
U = 2.28 nH, R l = 50 kfi (SQUID), C 2 = 0.5 pF, L 2 = 3.2 
nH, R 2 = 30 kfi, C 3 = 2 fF, L 3 = 10.8 nH, R 3 = 3 kfi, R = 
100 Q and T= 100 mK. 



spaced as would be the case for a classical inductance. 

In order to find a better quantitative agreement with 
the level spacing, the Schrodinger equation J2J was nu- 
merically solved under the assumption that the particle 
is localized in one of the wells. The experimental peaks 




FIG. 3: The position of the main resonances as a function of 
applied flux for the 4-SQUID sample and the calculated tran- 
sitions. The main resonances due to the SQUID-loops consist 
of a double-peak structure. The potential well and the energy 
levels and shown in the inset. The arrows indicate the transi- 
tions that are clearly seen in the experiment. The transition 
from first to fourth harmonic is shown by the dashed line. 

together with the calculated transitions are shown in Fig 
|3J The form of the cosine potential decreases the level 
spacing from the harmonic case. This deviation from the 



harmonic oscillator case is largest for the third transition 
as can be seen in Fig. [3| where also the transition from 
the first to fourth harmonic level is shown for comparison. 

The parameter for Ej used in the calculation is about 
twice the value given by the A-B relation (see Table 1). 
Typically, Ej is expected to be renormalized downwards 
due to the low impedance environment, but in our case it 
is renormalized upwards. Similar disagreements between 
the A-B value and have been reported before 0, Hif . 
Our model does not, however, explain the double peak 
structure (see Fig. [2J found in all the IV-curves. This 
peak splitting is fairly constant over the whole measure- 
ment range but the position of the double peak is dif- 
ferent for the three main transitions and, as can be seen 
in Fig. |21 the multi-photon excitations given by P{E)- 
theory are not (except for transition (1,1,0)) consistent 
with the observed peaks. 

The 2-point measurements with both one and two 
SQUID(s) showed similar behavior as the 4 SQUID mea- 
surements. The number of SQUIDs in the sample did 
not seem to have any significant effect on the IV-curve. 
Rather, there are notable differences between the 2-lead 
and 4-lead samples. In the 1-SQUID sample with only 
two leads the current dropped very fast when tuning 
down Ej (from 360 pA at $/$ = to 40 pA at 
$/$o = 0.4). This behavior can be explained when 
considering that the current through the circuit is given 
by two rates: the excitation of oscillator modes in the 
SQUID and their subsequent relaxation, which depends 
on the environment seen by the SQUID. In the 2-lead cir- 
cuits the current is limited by the down relaxation and 
the effect can be approximatively explained with the for- 
mula 



T i = 2^Re{Y(E ln // l )}(E ln /fi)R Q | (L>|n) | 2 



(3) 



where the admittance is given by Y(u) = [1 / '{iu)C 'net) + 
i?o] _1 and Coet is the capacitance of the detector junc- 
tion (0.8 fF) in series with the resistance of the environ- 
ment R (100 0). 

The position of the clearest flux-dependent peaks for 
the 1-SQUID sample are shown in Fig. 0] In this 
case, better agreement with the measured resonances is 
found when considering the full periodicity of the cosine- 
potential in Eq.|21 However, the value of Ej was taken to 
be twice the A-B value, as in the case with four SQUIDs. 
This apparent enhancement of Ej is probably due to the 
charging energy, as discussed in Ref. The effect ac- 
cording to the theory is, however, smaller than what we 
observe. 

Because the transitions are due to transfer of Cooper- 
pairs in the detector junction, the allowed first-order 
transitions can be found by calculating the matrix ele- 



ment | (I | e" 



lip 



between bands I and n. In addition, 



due to van-Hove like singularities, the observed transi- 
tions are between band edges. The transition between 
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the first and fourth band is clearly visible as two distinct 
peaks. The lowest bands are so narrow that they only 
show up in the width of the resonance peaks. As the 
Ej/Ec ratio is tuned down, the life-time of the states 
grows and this should lead to narrower peaks. But, in- 
stead we observe a broadening of the resonances, indicat- 
ing a broadening of the bands as expected from theory. 

The theory for Josephson junctions ^3] tells that in 
order for band formation we need to suppress the ohmic 
or quasiparticle dissipation in the environment, which 
causes the phase to localize. In our system, this sup- 
pression is provided by the large quasiparticle resistance 
of the detector junction. Therefore, the wavefunctions 
are 47r-periodic and each band from the 27r-periodic case 
is split into two. The observed transitions are, however, 
the same as what would be expected for 27r-period bands. 
Consequently, the true periodicity of the bands cannot be 
resolved in the experiment. 

In summary, we have experimentally studied the quan- 
tum mechanical energy levels of the Josephson junction. 
Our results for samples with Ej/Ec 3> 1 can quali- 
tatively be described by P(£')-theory, consistent with 
multi-photon excitations in the experiment. The non- 
linearity of the SQUID systems, prominent of Ej ~ Ec, 
can be taken into account by considering the exact form 
of the cosine-potential. Evidence of the existence of 
Bloch bands is observed in our 2-lead samples both in the 
form of van-Hoven like singularities between band edges 
and a broadening of resonance peaks. Our results show 
that a small superconducting junction can be employed 
as a detector for mesoscopic quantum circuits. 
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FIG. 4: Position of the measured resonances for the 1-SQUID 
sample and the theoretical transitions between band edges 
when the wavefunctions are 47r-periodic. The resonances are 
grouped into 4 groups, which are transitions between bands 
l-»2, l->3, l-> 4 and 1— 1-5. As in the 4-SQUID sample, 
there is a double peak structure which it not explained by the 
model. 



